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Abstract 

The interaction of an interfacial crack with small impurities is analysed on the basis of an asymptotic 
formula derived by the authors. The interaction between the main crack and the defects (e.g. small 
cracks or inclusions) is described asymptotically by analysing the dipole fields and the corresponding 
dipole matrices of the defects in question. The method is generic, and it serves interfacial cracks 
with general distributed loading on the crack faces, taking into account possible asymmetry in the 
boundary conditions, and in a particular configuration for a crack in a homogeneous medium results 
agree with those obtained earlier by [3|. Shielding and amplification effects of the defects on the 
propagation of the main crack along the interface are investigated. Numerical computations based 
on the explicit analytical formulae show potential applications in the design of composite and fiber 
reinforced materials. 
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1 Introduction 



Asymptotic models of a brittle crack interacting with a small defect (or a finite number of micro- 
defects) have been addressed in ([TJ [TSJ [T7J [TBI HI El HHJ HO]). A finite crack interacting with a large 
number of micro-defects in an elastic solid was also considered in ([7j[TT]); a survey of publications on 
macro-microcrack interaction problems is included in [13 . The approach [TTj has led to homogenization 
approximations for dilute composite media. Efficient numerical algorithms for cracks in solids have been 
developed for formulations based on singular integral equations. In particular, numerical techniques have 
been adopted to solve macrocrack-microdefects interaction problems f |15L \W\ 15]). The papers (|12L Ifj]) 
led to the notion of shielding/amplification due to the presence of micro-defects. A formal iterative 
procedure developed for an integral equation formulation describing interaction between a finite size 
interfacial crack and near-interface smaller cracks is discussed in [211 [55] ; the solutions of related integral 
equation involved the Chebyshev polynomial representations. Interaction of a dislocation dipole with a 
Mode III interfacial crack was studied in |2J; the results included the analysis of the "shielding effect" 
for the case when the dislocation dipole is placed in a neighbourhood of the crack tip. The classical 
mathematical model for bi-material media containing cracks propagating along the interface between 
different phases has been developed by |23| . 

In this paper we consider structures consisting of two different elastic materials, with the geometry 
and main notations shown in Fig. [TJ It is assumed that the middle point Y of the micro-defect is situated 
at a distance d from the tip of the main interfacial crack and the micro-defect length 21 is much smaller 
compared to d (that is. e = l/d is a small parameter, < e <C 1). 

Our asymptotic procedure can be summarised as follows. As e — > 0, the asymptotic solution is sought 
in the form 

u(x) = u<°> (x) + eW^ ( + e 2 « (2) (*) + O (e 3 ) , (1) 

where x — (xi^x^). The first term v^ ' corresponds to an unperturbed problem for a single interfacial 
crack. The second term denotes the boundary layer appearing near the small defect. The third 

term v& is the solution of the problem for the unperturbed solid with the interfacial crack loaded along 
the crack surfaces by tractions induced by the boundary layer 

For a particular case of a crack in a uniform elastic body and the remote load applied away from the 
crack tip, our asymptotic approximation agrees with that of |3J: 

cos — a J cos I — — a J + O (e 1 . (2) 



K m 4 V 2 / V 2 

The procedure outlined above is the essence of the so-called general dipole matrix approach developed 
in (TJJ and described for various types of small defects in Here, we employ the representation for an 
interfacial crack developed in [T3] where the subscripts ± correspond to the position of the point under 
consideration within dissimilar half-planes (see Fig. [TJ. 

The present paper addresses the accuracy of the proposed asymptotic approximation and it analyses 
the discrepancy in the traction transmission conditions on the interface due to the boundary layer W^ 1 
and hence its contribution to the perturbation of the stress intensity factor. 

2 Preliminary results (unperturbed crack) 

Using the weight function approach, the stress intensity factor for the interfacial crack loaded by given 
out-of-plane tractions p^(xi) (x\ < 0) can be evaluated as 

4°) = fl[° »-P + ^+^P~^\ - Xl) -^ dxv ( 3 ) 
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Figure 1: A line defect (microcrack/rigid line inclusion) of length 21 arbitrarily oriented and situated at a distance d 
from the tip of the main interfacial crack. Fig. Qa) corresponds to the problem under consideration. Fig. Qb) illustrates 
a configuration where the remote load is characterised via the value of the stress intensity factor. 



Without loss of generality, we formulate the results in terms of concentrated forces applied on the 
crack surfaces, so that the formula ^ can be rewritten in the following form: 



(0) 

in 



1+7? 

V2n 



-i/2 1 — rj 



-/ -\-l/2 



(4) 



Here r\ is the contrast parameter, r\ = — + jtz+), and the point forces F^ are applied at the 

points (— a~j,0+), (j = 1, ...,n), while F~ are applied at the points (— aj ,0— ), (j = 1, ...,m). 
The formula Q can then be written in the form 



k (0) - V x (0 ' +) + Vir (0 ' _: 



(5) 
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We note that the stress intensity factor produced by a point force F 1 ^ applied to the upper /lower 
crack surface at a distance ar behind the crack tip has the form: 



A', 



(o,±) _ l±»7 p± / ± v-l/2 



(0) 



Correspondingly, the components of the displacement gradient computed at an arbitrary point Y± 
(dcos0,dsin^>) can be evaluated for the point force Fj applied at the point (— a+,0+) as follows 



3u(o,i) 
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(7) 



0u(Oj) 



dx 2 



i 7 



nd(fi + + /j_)(2cos0 + a^/rf + d/aj) 



1 / a 

2 h 



d 



, ■ cos — 
d 2 



d 3(j) 
■cos — 
a, 2 



(8) 



4 



and for the point force F- applied at (— a ■ , 0—) we have 



Y± nd((x + + fi-)(2cos(f> + cij /d + d/dj ) 
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Here the superscript ± corresponds to the surface where the point force F, is applied. 



(10) 



3 Refined asymptotic formulae 

The perturbation of SIF produced by a linear defect (microcrack/rigid line inclusion) of length 21, centred 
at the point Y = (d cos 0, d sin 0) , making an angle a with the positive £i -direction and situated away 
from the interfacial crack tip (l/d <C 1) was obtained by [9]: 



AK 



in 



2 /X+M " B(d, 0) • M(l, a)c(d, 0) + O ( ~ 



where 



B(d,0) = £ 
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and A4(Z, a) is the dipole matrix: 
- for micro-crack 



- for rigid line inclusion 



M cr (l, a) — —irl 2 
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(11) 
(12) 

(13) 
(14) 



Finally, in the presence of several micro-cracks/inclusions, formula (111 will contain a sum over all 
defects. 



Formula (111 can be written in the form similar to that of (|2j 

where the function G is different for the micro-crack and the rigid line inclusion: 

y/nd 



,30 
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(16) 
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G. m = sin (~ - a] nd B(d, 4>) ■ [cos a, sin a]. (17) 
V 2 J y/2K$l 

We would like to note that in the case when a — 30/2 = ±7r/2 respective formulae Q as well as ( 15 1 and 
(16) provide the same (zero) results regardless of particular values of all other parameters. The same is 
valid for the rigid line inclusion when 30/2 — a = or ±7r. 

In what follows, we compare the numerical findings and asymptotic approximations and comment 
on applicability of the results of [3] . 

First, we compare numerically formula ([2| with (15 1 in the case when the main crack lies in a 
homogeneous plane (fi + = /x_ = //) and a microcrack of the length I = 0.01 is situated at the distance 
d = 1 from the crack tip (that is, e = 0.01). A symmetrical loading consisting of two point forces 
F + = F~ = F acting at a distance a from the crack tip is applied. 

In Fig. [2] we present the relative error for the case when one applies ^ instead of ( Jl5[ ). Other param- 
eters involved in the computations are: a — 0; 7r/3.9; ir/2 (cases 1,2,3 respectively), <fi = 7r/4, ir/2, 37r/4 
(cases a, b, c respectively). 




Figure 2: Relative error of the formula (J2J) in comparison with more accurate formula |l5| as a function of the dimension- 
less parameter a /din the case of a homogeneous plane with the main crack and a micro-crack situated at a distance d from 
the crack tip. A symmetrical point loading is prescribed on the crack surfaces at a distance a. Other parameters a and <f> 
describe the position and orientation of the micro-crack (see Fig.[TJ. The following values were used in the computations: 
a = 0; 7r/3.9; it/2 (cases 1,2,3), <j> = 7r/4, w/2, 37r/4 (cases a, b, c respectively). 



One can see there is a pronounced contribution from the load applied close to the crack tip. Only if 
the ratio a/d is big enough the asymptotic approximation Q is sufficiently accurate (for example, for 
a/d > 100 the relative error is less than 1%). 

For a separate point force applied on the upper (F + ) or lower (F~) crack surface we have the following 

asymptotic formula (with the accuracy O (^(d/a) 3 ^ 2 
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(18) 



Using these representations we can perform asymptotic analysis for an arbitrary number of points 
forces to obtain the final result for SIFs in the case of a single defect situated at the point Y±: 
- for a micro-crack 
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for a rigid line inclusion 
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(19) 



(20) 



The terms XI1X2 can be estimated as O (max{d, b}/a). This result shows that the second term of the 
order 0(Z 4 /d 4 ) in the formula Q is significant if the following estimate is valid max{d, b}/a = o(l 2 /d 2 ), 
otherwise there are other essential terms in the asymptotics larger than 0(l 4 /d 4 ). 



Note that, when computing (19 1 and (20 1, the first terms in (18) have disappeared after summation 
as the load is self-balanced. 



To give an additional comparison with (19 1, we also consider a special case when the crack faces are 
loaded by a "three-point" loading system consisting of a point force F acting on the upper crack face at 
a distance a behind the crack tip and two point forces F/2 acting on the lower crack face at distances 
a — b and a + b behind the crack tip. The parameter b is fixed. We assume here that b/a = o(d/a) and 
obtain 
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for a microcrack, while for a rigid line inclusion it takes the form: 
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4 Numerical computations 

Before the discussion of computations we would like to point out that the simplified asymptotic formulae 



(19) and (20 1 give similar predictions regardless of the bimaterial parameter 7/ (the ratio of the shear 
moduli /i + //i_). Of course, the value of the stress intensity factors will be a function of this ratio. The 



more accurate formula (15 1 indicates the influence of possible non-symmetry of the applied load. 

We first present the results for a homogeneous body (/i + = /i_ = /i), see Fig. [3} In the computations, 
we use the value of the small parameter e = 0.01 which describes a relative size of the small defect in 
comparison with the distance to the crack tip (d = 1 and / = 0.01). 

On the diagrams, the horizontal axis stands for the angle cf> {<ft £ (— 7r, 7r), see Fig. [T| defining the 
position of the center of the micro-defect with respect to the crack tip. On the vertical axis we measure 
the value of the angle a (0 < a < n). The diagrams show the borders between the region where AK < 
(shielding effect) and AK > (amplification of the SIF due to the presence of the small defect). The 
respective regions are shadowed by light gray (shielding) and dark gray (amplification). 

Fig.|3ja)-upper corresponds to a two-point (symmetrical) load situated at a distance a from the crack 
tip in the case of a small micro-crack placed at the distance d from the crack tip. As it is clear from 
the presented results, the distance between the load support and the crack tip plays an important role 



7 



(a) micro-crack 



(b) rigid line inclusion 



Figure 3: " Shielding" (light grey) and "amplification" (dark grey) regions created by a micro-crack (a) or a rigid line 
inclusion (b), whose position relative to the xi-axis is characterised by the angles <j> and a. Upper part: two-point loading 
(ft = 0) at a distance a from the crack tip. Lower part: non-symmetrical three-point loading (ft = 1) at a distance a from 
the crack tip (see Fig. Qa)). The borders between the regions for a = 2, a = 4 and Gong's approximation are shown by 
solid, dotted and dashed lines, respectively. 

(solid, dotted and dashed lines are different in the middle of the diagram). Moreover, in the case a = 2 
when the load is close enough to the crack tip (but still a > d) new small amplification regions appear 
in the corners of the picture. It is natural that the diagram is completely symmetrical with respect to 
the angle (f> as the load is symmetrical and the material parameters do not influence the diagram. 

The results change when we consider a three-point non-symmetrical load (as shown in Fig.JTJa)). In 
Fig. [3ja)-lower we set the parameter b, the distance between two point forces on the lower crack surface, 
to 1. 



In Fig. [3jb)-upper we consider the same symmetrical load as in Fig. [3]ja)-upper, however now a small 

defect is represented by the rigid line inclusion placed at the same distance d from the crack tip. The 
shielding and amplification regions are completely different in comparison with the micro-crack case. As 
in Fig. [3ja)-upper, due to the symmetry of the load, the diagram of Fig. [3^b)-upper is symmetric with 
respect to the angle 4>. 



In Fig.|3jb)-lower we consider a three-point non-symmetrical load as in Fig.j3j[a)-lower, but we replace 

the micro-crack by a rigid line inclusion. As in the previous cases, when a is large (and even for a = 4) 
non-symmetry in the load is hard to observe. However, in the case a = 2 the diagram shows a clear 
asymmetry. 

Finally we would like to discuss the influence of the inhomogeneity on the shielding effect. As we just 
mentioned, it can be only observable for a highly pronounced non-symmetrical loading which is applied 
close to the crack tip. For this reason in Fig. [4] we consider the aforementioned three-point loading for 
the case when b = 1.9 and a = 2,4 and 10000. Two ratios of the material parameters were used in the 
computations: = 0.1 (upper part) and = 10 (lower part). 



The diagrams of Fig. Qa) correspond to a micro-crack defect, while those of Fig. [4|b) correspond to 

a rigid line inclusion. 






(a) micro-crack 



(b) rigid line inclusion 



Figure 4: "Shielding" (light grey) and "amplification" (dark grey) regions created by a micro-crack (a) or a rigid line 
inclusion (b), whose position is characterised by the angles <f> and a (see Fig. [TJ. The diagrams are for a three-point loading 
with b = 1.9, and the material parameters ratio = 0.1 (upper part) and fi-^/fi- = 10 (lower part). The borders 

between the regions for a = 2, a = 4 and a = 10000 are shown by solid, dotted and dashed lines, respectively. 



5 Discussion and conclusions 

In this paper we have discussed asymptotic formulae for the perturbation of the Mode III SIF due to 
the presence of a small defect at a large distance from the tip of an interfacial crack (e = 1/d « 1). 
We have derived a novel accurate analytical formula based on the dipole matrix and weight-function 
approach. It was noted that the position of a small defect influences the SIF leading to shielding and 
amplification effects for various geometrical and material parameters. We have also shown that, for a 
Mode III deformation, the material inhomogeneity (dissimilar plane) plays an important role only in the 
case of a non-symmetric loading applied close to the crack tip. 

The method exploited here allows to define Mode III SIF for an arbitrary load by generalising the 
asymptotic representation and to estimate the error of the asymptotic approximation. Finally, we 
would like to emphasise that, due to the linearity of the problem, the influence of a finite number 



of small defects can be computed by a simple summation of the formulae (15 1 over every particular 
defect. The dipole matrices for clusters of small defects would have to be employed if the number of 
micro-cracks/inclusions becomes large. 

Compared to the earlier publications on related topics (e.g. [T3], [2T], [55], 0), our results allow for 
an explicit analytical representation of the asymptotic solution. Also, the approach based on the analysis 
of dipole fields is generic, and it is not restricted to a particular shape of defects, like micro-cracks, for 
example. The method is applicable equally well to the case of small voids or inclusions of arbitrary 
shapes, with smooth or non-smooth boundaries, including clouds of defects of different shapes. 
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